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Abstract. We study the quark-hadron mixed phase in proto-neutron stars with the finite-size 
effects. In the calculations of pasta structures appeared in the mixed phase, the Gibbs conditions 
require the pressure balance and chemical equilibrium between two phases besides the thermal 
equilibrium. We find that the region of the mixed phase is limited due to thermal instability. 
Moreover, we study the effects of neutrinos to the pasta structures. As a result, we find that the 
existence of neutrinos make the pasta structures unstable, too. These characteristic features of 
the hadron-quark mixed phase should be important for the middle stage of the evolutions of 
proto-neutron stars. 



1. Introduction 

There are many uncertainties left for the hadron-quark phase transition, e.g. the equation of 
state (EOS) of quark matter or deconfinement mechanism. Assuming the quark deconfinement 
transition to be of the first order, since there are two conserved quantities, baryon number 
and electric charge, the phase equilibrium in the mixed phase must be carefully treated by 
applying the Gibbs conditions [1], instead of the Maxwell construction. A simple treatment of 
the mixed phase may be the bulk Gibbs calculation, where phase equilibrium of two bulk matter 
is considered without electromagnetic interaction and surface tension. Generally the properties 
of the mixed phase should strongly depend on electromagnetic interaction and surface tension, 
and these effects, sometimes called "the finite-size effects", lead to the non-uniform "Pasta" 
structures. The EOS of the mixed phase becomes similar to the one by the bulk Gibbs calculation 
for weak surface tension, and to the one given by the Maxwell construction for strong surface 
tension [21 El SI IS] - The charge screening is also important for their mechanical instability. In this 
paper, we study the hadron-quark mixed phase with the finite-size effects at finite temperature 
by extending the previous works [H El [6] . We also present some results for the neutrino-trapping 
case and discuss its effect. 



2. Equation of state for each phase 

The theoretical framework for the hadronic phase of matter is the nonrelativistic Brueckner- 
Hartree-Fock approach including hyperons such as £~ and A [7J. There is a controversy about 



the Y<~-N interaction. The recent experimental result about hypernuclei has suggested that it 
is repulsive [HJ|9], while we use here a weak but attractive interaction; S~ then appears before 
A in uniform and neutral matter. It would be interesting to see how our results are changed 
by using different Y,~-N interaction, and we will discuss it in the future work. Moreover, we 
adopt the frozen-correlation approximation at finite temperature [10] . This approximation is 
feasible for temperature T < 50 MeV. 

For the quark phase, we adopt the MIT bag model or the density dependent bag model 
consisting of it, d ,s quarks [11]. Probably, these models are too simple to describe quark 
matter in a realistic way. We will adopt more sophisticated models in the future [12] . We 
assume massless u and d quarks and s quarks with the current mass of m s = 150 MeV. We set 
the bag constant B to be 100 MeV fm~ 3 in the MIT bag model. In the density dependent bag 
model, we set the vacuum energy B(p) as 

B( P ) = £00 + (B - Boo) exp[-/3-] 

Po 

with Boo = 50 MeV fm -3 , Bq = 400 MeV fm -3 , and j3 = 0.17. For a more extensive discussion 
of this topic, the reader is referred to Nicotra et al. [11] and references cited therein. 

For the mixed phase, we assume non-uniform matter with various geometrical structures 
("pastas"); droplet, rod, slab, tube, and bubble. Introducing the Wigner-Seitz cell to treat such 
non-uniform structure, we use the local density approximation for particles. Between the hadron 
and quark phases we put a sharp boundary with a constant surface tension parameter. We 
impose the Gibbs conditions, i.e. chemical equilibrium among particles in two phases consistent 
with the Coulomb potential, and a pressure balance consistent with the surface tension. Here, 
we write down the conditions of the chemical equilibrium in the mixed phase; 

Pu + Pe - Pv £ = fJ-d = Ps, 

fip + fi e - fl Ue = [i n = HA = fi u + 2fi d , (1) 
MS" + Pp = 2 Pn- 

Note that we do not take into account the anti-particles in this paper. 

Since the role of the surface tension to the mixed phase has been already studied in the 
previous papers [U [5], we use one constant value of the surface tension parameter as = 40 
MeV fm -2 in this paper. 

Finally, we compare the Helmholtz free energy among the pasta structures and choose which 
structure is most favored. The Helmholtz free energy for each cell is then given as 

F = f dr 3 F H M + f dr 3 F Q [n q ] + F e + F Ue + E c + E s (2) 
Jv„ Jv Q 

with i = n,p, A, a, q = u, d, s, Fh (Fq) is the Helmholtz free energy density for hadron (quark) 
matter, and Es = crS the surface energy with S being the hadron-quark interface area. F e and 
F„ e are the free energies of the electrons and neutrinos. For simplicity, muons are not included 
in this paper. 

3. Results 

In proto-neutron stars, we must consider mainly two effects, namely temperature and neutrino 
fractions. The range of temperature is roughly ~ - 40 MeV, and the range of neutrino fraction 
is ~ - 0.2 [13]. In previous paper, we have studied the effects of temperature on the quark- 
hadron mixed phase [6]. In this study, we check the effects of neutrino fractions and summarize 
the finite size effects of the quark-hadron mixed phase in the proto-neutron stars. 



3.1. Effects of temperature 

In previous study, we have shown that finite temperature makes the mixed phase unstable [6]. 
FigureUJis one example to show the mechanical instability of the pasta structures as temperature 
is increased. Here, we assume that the vacuum energy is constant, B = 100 MeV fm -3 , applying 
the MIT bag model. The left panel of Fig. [1] shows the free energies per baryon of the droplet 
structure for several values of temperature. The quark volume fraction in each cell (R/R\y) S 
is fixed to exclude the trivial R dependence. Here we use the optimal value of (R/Rw) 3 at 
T = MeV in every curve. We normalize them by subtracting the free-energy at infinite radius, 
AF = F(R) — F(R — > oo), to show the R dependence clearly. The structure of the mixed phase 
is mechanically stable below T ~ 60 MeV, but the optimal value of the radius R is shifted to 
the larger value as T increases. This behavior shows a signal of the mechanical instability and 
it comes from the charge screening effect and the thermal effect. To elucidate this point more 
clearly, we show each contribution to F/A in the right panel of Fig. [TJ i.e., the Coulomb energy 
per baryon Eq/A, the surface energy per baryon Eg/ A, and the correlation energy per baryon 
E corr /A. The baryon density and the temperature are set as ns = 2 no and T = 50 MeV. The 
correlation energy is produced by the difference of the bulk energy due to the rearrangement of 
charged particles. The entropy term in the bulk energy becomes dominant over the Coulomb 
energy in this case to exhibit the instability. The details about the effects of temperature, 
including EOS, are fully discussed in our previous paper [6]. 




Figure 1. The left panel shows that the droplet radius R dependence of the free-energy per 
baryon for n# = 2 uq and different temperatures. Here, we do not take into account neutrinos 
to see the effects of temperature. The quark volume fraction (R/R\y) 3 is fixed to be the optimal 
value at T = MeV for each curve. The free energy is normalized by the value at R — > oo. The 
filled circles on each curve shows the energy minimum. The results are for B = 100 MeV fm -3 , 
a = 40 MeV fm -2 . The right panel shows each contribution to the R dependence of the free 
energy, the Coulomb energy, the surface energy, or the correlation energy per baryon (E c /A, 
E 8 / A, E corr /A) at T = MeV (thick lines) and T = 50 MeV (thin lines). 



3. 2. Effects of neutrinos 

The left panel of Fig. [2] shows the free energies per baryon of the droplet structure at several 
values of neutrino fractions. Here we use the optimal value of (R/R\y) 3 at Y Ue = in every curve. 
In these figures, we adopt the density dependent bag model shown as Eq. (JJ). We normalize 
them by subtracting the free-energy at infinite radius, as same as Fig. [TJ The structure of the 
mixed phase is mechanically stable below Y Ve ~ 0.10, but the optimal value of the radius R is 
shifted to the larger value as Y Ue increases. 



To elucidate this point more clearly, we check each contribution to F/A in the right panel of 
Fig.[2]again. The baryon density and the temperature are set to be ub = 2.5 no and T = 10 MeV. 
In this, we compare the case for Y Ve = 0.01 (thin lines) with the one for Y Ve = 0.15 (thick lines). 
We can see the main contributions responsible to the change of free energy are the correlation 
energy and the Coulomb energy, though the thermal effect appears only in the correlation energy 
in Fig. [TJ 

Why does the Coulomb energy become important for the neutrino trapped case? With the 
presence of neutrinos, the electron fractions becomes rich because of the chemical equilibrium 
Eq. @. Accordingly the fraction of particles with positive charge is enhanced through the 
charge neutrality condition. Thus, the net charge density is reduced over the whole region of 
the cell. Consequently the Coulomb energy is largely reduced with neutrinos. The change of 
the correlation energy follows that of the Coulomb energy since the rearrangement effect is also 
reduced. Both the Coulomb and correlation energies give rise to the instability in this case. 





Figure 2. Same as Fig. [H but for different neutrino fractions Y Ve at constant temperature, 
T = 10 MeV. In this case, we set baryon density as ub = 2.5 no, and the quark volume fraction 
(R/Ry/Y is fixed to be the optimal value at Y Ve = 0.01 for each curve. The thick lines of the 
right panel are the case of Y Ve = 0.15, the thin lines Y Ue = 0.01. 

Fig. [3] shows the change of the Coulomb energy by trapped neutrinos more clearly. The figure 
shows the density profiles within the ID cell (slab) for ng = 2.5 no at Y Ue = 0.01 (left panel) and 
Y Ue = 0.01 (right panel). Each temperature is set to be T = 10 MeV. The volume ratio (R/Ryy) 3 
of the quark phase is fixed to be the optimal value at Y Ue = 0.01. Clearly, the magnitude of 
the Coulomb energy is decreased at high neutrino fraction. At high neutrino fraction, since 
the enhancement of electrons increases the particle fractions with positive charge, the negative 
charge of the quark slab is suppressed. As the result, the magnitude of the Coulomb energy by 
the quark slab is decreased. Of course, the particle fractions with positive charge increase in the 
hadron slab, too. But, they do not change so much because of their heavy masses. 



4. Summary and Discussion 

We have studied the hadron-quark phase transition and the properties of the mixed phase in 
proto-neutron stars. We have taken into account the finite-size effects by imposing the Gibbs 
conditions on the phase equilibrium, and calculating the density profiles in a self-consistent 
manner. 

We find that the hadron-quark mixed phase become unstable at high temperature and/or 
high lepton fraction. The instability by finite temperature mainly comes from the correlation 
energy, while the effects of neutrinos on free energy appear both in the correlation energy and 
the Coulomb energy. 
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Figure 3. Density profiles and Coulomb potential (j> within ID (slab) for ub = 2.5 tlq at T = 10 
MeV. Here, the neutrino fractions are set to Y Ue = 0.01 (left panel) and Y Ue = 0.15 (right panel). 
The cell sizes are R\y = 30 fm in these figures. The slab size are R = 10.9 fm. (R/Rw) 3 is fixed 
to be the optimal value at Y Ue = 0.01. 



Hence, we can suggest that the quark-phase transition is close to the EOS with the Maxwell 
condition at the very first stage of proto-neutron stars, because the temperature is high, T ~ 40 
MeV, and the neutrino fractions are also high, Y Ue ~ 0.2. However, in the middle stage of 
evolution of proto-neutron stars, the pasta structure will appear, where initial cooling and 
deleptonization proceed. Hence, the simple Maxwell condition might be enough to take into 
account the quark-hadron phase transitions in the simulations of supernovae or black hole 
formations, although such treatment is inappropriate in the simulations of proto-neutron star 
evolutions. We can see the similar discussions in the recent papers by Schaffher's Group, though 
they have not included the finite-size effects [21 [15] . 

Finally, we note again that the EOS has many uncertainties, especially for quark matter. We 
simply adopted the MIT bag model and the density-dependent bag model in this paper, while 
other quark models, including chiral restoration or color superconductivities, may change our 
results [12]. These are open questions for astrophysics and nuclear physics. 
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